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The gas-liquid transition is a first-order transition terminating at a finite-temperature crit- 
ical point with diverging density fluctuations. Mott transition, a metal-insulator transition 
driven by Coulomb repulsion between electrons, has been identified with this textbook tran- 
sition. However, the critical temperature of the Mott transition can be suppressed, resulting 
in unusual quantum criticality. It accounts for non-Fermi-liquid-like properties, and strongly 
momentum dependent quasiparticles as in many materials near the Mott insulator. Among 
all, the mode-coupling theory of the density fluctuations supports d-wave superconductivity 
at the order of a hundred Kelvin for relevant parameters of the copper oxide superconductors. 
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Metal-insulator transitions are characterized by a huge change of electrical conductivity. 
Among all, the Mott transition is driven by the electron Coulomb repulsion. 1,2 By using the 
dynamical mean-field theory, 3 Kotliar et al. suggested that this transition is equivalent to the 
gas-liquid transition and thus to the textbook transition of the Ising model. 4 ' 5 Pressure control 
experiment on vanadium oxides supported the equivalence. 6 As in the gas-liquid transition, 
the electron density susceptibility diverges at the critical point. 

Recent experimental results near the Mott insulator have clarified diversity such as renor- 
malized Fermi liquid, charge and magnetic orders, and superconductivity with their compe- 
titions. 2 In spectroscopic data, strongly momentum-dependent structure appears as another 
conspicuous feature particularly in the copper oxide superconductors. 7,8 A well known exam- 
ple is the flat quasiparticle dispersions in the (n, 0) and (0, 7r) regions of the cuprates clarified 
in the photoemission experiments. Spatial inhomogeneity of electrons is also suggested in the 
cuprates 9, 10 and the manganites. 11-13 Severe competitions and self-organized heterogeneous 
structure in real and momentum spaces are a way of characterizing the region near the Mott 
insulator. The relation of these remarkable features to the above mentioned Mott transition 
but crucially modified by quantum effects was discussed recently. 14 

In this letter we first extend the previous considerations 14 especially in two dimen- 
sions. 3,6, 16,17 The critical exponents of the Mott transition thus obtained indeed agree with 
the recent experimental results on a K-(ET)2-type organic compound. 15 As a main result, we 
show that the mode coupling theory supports the d-wave superconducting state driven by this 
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quantum Mott criticality, at the order of a hundred Kelvin for realistic parameter values of 
the copper oxide superconductors. 

In the Ginzburg-Landau (GL) theory 18 of the gas-liquid transition, the free energy is 
expressed as F = A(T)n 2 /2 + Bn /4 with A oc T — T c , and a positive constant B. The natural 
order parameter to be controlled is the relative density n measured from the critical point. At 
the critical point with the temperature T = T c , the density susceptibility (or compressibility) 
defined by \ c = [d 2 F/dn 2 \~ l diverges. Below T c , the spinodal decomposition or nucleation 
drive the system into a simple equilibrium with phase separation or coexistence. 

In the Mott transition, one can control either the electron filling n (equivalently the hole 
density X = 1 — n) or the bandwidth (or transfer) t relative to the local Coulomb repulsion 
between electrons, U. When U/t increases, the conjugate variable to U, namely the doublon 
density D signals the transition by discontinuous or singular reduction in £). 3,19 ' 20 Here, a 
doublon is defined by a site with a spin-up and a spin-down carrier electrons coexisting. In 
terms of the natural order parameters, £ = X or D, the transition has a similar anomaly to 
that of the gas-liquid transition. 4 ' 5 We employ the notation of the conjugate variables to D 
and X as [id = U and the chemical potential fix = M> respectively. Although the full free 
energy may be expressed by a functional of the Green function, 21 the essence is expressed 
by these natural order parameters after integrating out other degrees of freedom. At T c , 
exactly as in the gas-liquid transition, the density susceptibility (or charge compressibility) 
Xx = [d 2 F / dX 2 ]^ 1 and the doublon-density susceptibility xd = [d 2 F/dD 2 ]~ l diverge, below 
which the first-order transition line with a jump of £ = X or D starts. Since X is conserved, 
the filling-controlled first-order transition appears as the phase separation. 

Additional factor for realistic electrons is the charge neutrality. Even when T is below 
T c in the Hubbard-type models considering only the short-ranged force, in reality, the phase 
separation dynamics freezes at the stage allowed from the long-ranged repulsion, although the 
homogeneous phase is still unstable in the spinodal region. This stabilizes a spatially inhomo- 
geneous structure with the length scale determined from the balance of the spinodal instability 
to the electrostatic condition. 14 At T > T c , if impurity potential or lattice distortions exist, 
it may also drive the inhomogeneity because of the enhanced density susceptibility. This re- 
minds us of the structure observed by scanning tunnel microscope (STM) in the cuprates and 
manganites. 9-13 The present results obtained from the Mott criticality has a tight connection 
to the approach from dynamical stripe fluctuations. 22 ' 23 

How do quantum effects modify this standard classical picture? In the usual quantum 
critical phenomena, with increasing quantum fluctuations g, T c becomes lowered to zero be- 
yond which the transition disappears as in Fig. la. On the contrary, the boundary between a 
metal and an insulator does not terminate at T = 0. At T = 0, metals and insulators in clean 
systems cannot be adiabatically connected and clearly distinguished by zero or nonzero Drude 
weight and density susceptibility. 14 ' 24 ' 25 Therefore we have the feature in Fig. lb, in which, at 
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Fig. 1. Comparison of conventional quantum transition (a) with the quantum Mott transition (b) in 
the parameter space of fi^, T and the parameter g to control quantum fluctuations. In the Mott 
transition (b), a zero-temperature critical line continues beyond the marginal quantum critical 
point (solid circle) even after the first-order transition surface (shaded sheet) shrinks. 

T = 0, the first-order and continuous transition lines meet at the marginal quantum critical 
point. 

We now further study implications of the marginal quantum criticality. A simple way of 
understanding the quantum effect is to start with the quasiparticle representation at T = 0. 
Except one-dimensional systems, the Fermi liquid theory describes metals by quasiparticles, 
where the interaction is renormalized to the single-particle coefficient. In the insulating side, 
the single-particle Green function G may also be given from a quasiparticle description with 
a nonzero gap A c as G(q, w)^ 1 = (u; + A c + E(q, k)) in the hole-doped region (, namely, for the 
pole of the lower Hubbard branch). Here, E(q,k) represents the expansion around the Mott 
gap edge in terms of k as E(q, k) = a(q)k 2 + b(q)k 4 + • • • , where k is the momentum perpen- 
dicular to the locus of ReE(q, k) = and q denotes that parallel to the locus KeE(q, k) = 0. 
Here Re E(q, k) > should be satisfied and vanishes at the gap edge. At the transition point, 
A c vanishes. The imaginary part of the self-energy, not considered here, does not alter the 
results on singularities below. If the dispersion at the Mott gap edge follows E oc k Zs , the 
dominant singular part of the free energy in terms of X is given by F oc x <yd+Zs ^/ d and hence 
the susceptibility xx oc X l ~ Zs l d for the spatial dimension d. 14 > 25 A similar scaling holds when 
the bandwidth is controlled by taking D measured from the Mott insulator value. Then, for 
a > at d > 2, the susceptibility does not diverge, while it diverges if a = at d < 3. Along 
T c = line, we expect a > 0, if xc i s finite. In contrast, at the first-order transition as well 
as at the marginal critical point, the susceptibility X( diverges. Therefore, when the marginal 
quantum critical point is approached at T = from the T c = line, the dispersion around 
the Fermi level has to become flat with a(q) = at least in a part of the Brillouin zone. This 
alters the dynamical exponent of the single-particle excitation, z s from 2 to 4, because the 
a term vanishes leaving quartic b term as the dominant one. The value a in reality depends 
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on momentum along the locus E(q, k) = 0. Near the marginal critical point, generically a 
at the point with the smallest amplitude becomes zero first. Then a quartic dispersion dom- 
inates at these special points, which results in diverging X(- Therefore, a singular electron 
differentiation generically evolves on the large Fermi surface, which consistently justifies the 
above scaling form of the free energy. This differentiation of electrons naturally explains a 
strong momentum dependence of quasiparticle dispersion observed in photoemission experi- 
ments with flat dispersions. 7 ' 8 The exponent z s = 4 and xx oc X^ 1 were supported in several 
numerical calculations in two-dimensional models when filling is controlled at T = 0, which 
suggests that the marginal critical point is close in these cases and a wide scaling region up 
toX~0.3. 2 < 27 > 28 

Now by rescaling numerical factors, the free energy at T = is expressed as F = — /i^C + 
a ^d+2)/d + h( -(d+A)/d + c( »(d+6)/d . . . with the constra i nt £ > 0. The metal-insulator transition 

across the T c = line with a > is driven by the ^ term. The marginal critical point 
at T = may be reached at a control parameter g = g c , with a = ao(g — g c (T)), ao > 
and b > 0. Although we consider the "mean-field" behaviour, the d-dependent F at T = is 
different from the classical GL form and gives unusually d-dependent exponents. The exponent 
(3 defined by the order parameter at g < g c and fi^ = as ( oc \g — g c \@ is given by f3 = d/2. 
Near g = g c , X( is expressed as xc = [d 2 F/dC, 2 Y l ~ [2(d + 2)a (g - g c )( 2,d ~ 1 /d 2 + 4(d + 
4)&C 4/d_1 Id 2 + 6(d + 6)cC 6 / d-1 /d 2 ] -1 - Then X( oc {g - g c )^ holds in the metallic side g < g c , 
yielding the "mean-field" exponent 7 = 2 — d/2. At g = g c , xc ~ & _1 C 1_<5 with S = 4/d is 
obtained. The scaling relation (55 = 7 + f3 is satisfied. We note that the scaling relation and 
the exponents are also derived from the free energy F(a,fi^) = S,~ d ~ Zt f(a^ y9 , fi^ yfl ) with a 
scaling function / and the correlation length £ oc (g — g c )- l l 2 oc which means y g = z. 

Here the crossover exponent = 4 is derived from the dynamical exponent of the density 
(two-particle) fluctuations given by zt = 4. When T c is low enough but nonzero, one may 
replace g — g c with T — T c in the above formalism, since T may also play a role of controlling 
quasiparticle dispersions. In two dimensions, it is reduced to F = a>o(T — T^g))^ 2 + 6£ 3 + c£ 4 
with T c (g c ) = and 

X C = [d 2 F/dC 2 Y l - [2ao(T - T c ) + 6b( + 12c( 2 ]- 1 . (1) 

Then, 7 = 1,(3 = 1 and 5 = 2 hold. 15 Remarkably, this agrees with recent experimental re- 
sult by Kagawa et al. 15 Very close to the critical point, the present mean-field exponent only 
marginally breaks down, because the Ginzburg criterion d+zt > (2(3+j)/u with v = 1/2 being 
the correlation length exponent indicates that the system is always at the upper critical dimen- 
sion. Here we restrict ourselves to the mean field study. Although the quasi-particle picture 
does not hold, the one-dimensional Hubbard model with next-nearest-neighbour transfers also 
satisfies the above scaling at its marginal critical point, giving 5 = 4. 26 At T > 0, we have the 
entropy term T( In ( in addition to the above F. This generates an essentially singular contri- 
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bution with the extremum at Co ex P[^c/^1) which does not contribute to the present scaling 
behavior near T = 0. Whereas at high temperatures without the Fermi degeneracy, the expan- 
sion around the extremum reproduces the regular GL form Ft = Aq(T— T C )Q 2 / A— fj,^, 
by redefining £ — ("o — ► C- 

Near the critical point, the density susceptibility is enhanced. We show consequences of 
the marginal quantum criticality within the mean-field scheme by assuming the dominant 
part of the susceptibility around a small and nonzero momentum Q as 

X ^ )= - l u J + D s {K^+( q -QY + ...)■ (2) 
with a diffusion constant D s , the distance from the criticality K and a constant V. They 
are scaled as T" 1 oc £ d ~ 4 ,K oc and D s oc £~ 2 , which reproduces the above scaling 
with z t = 4. By using this phenomenological form, Moriya has developed the self-consistent 
renormalzation (SCR) theory for spin fluctuations. 29 Now we formulate a similar scheme 
by considering the mode coupling for the electron-density and doublon fluctuations coming 
from the Mott criticality. Although we combine a perturbation treatment below, the essential 
ingredient of the Mott criticality may be incorporated by this scheme. We summarize the 
mode coupling method for density fluctuations. The inverse susceptibility at zero frequency 
is written from eq.(l) as 

xdQ, 0)^ = 2ao(T - T c ) + 6b(( Q ) + 12c^(C g 2 ) (3) 

i 

with the density fluctuation average 

1 f°° 

(C ' )= 2W ^ cot M^/2T)Im^(g,cu). (4) 

Eqs. (2)-(4) constitute self-consistent equations to account for the Gaussian fluctuations 
through the mode coupling. 29 The electron self-energy is calculated perturbatively for the 
filling control as 

£(<?, u n ) = ^2 G ( k > iu n )xx{q ~ k, i(u n - u m )). (5) 

k,n 

In the actual calculations, though some uncertainty remains, we choose the parameter values 
of a two-dimensional system appropriate for the copper oxide superconductors inferred from 
the frequency dependence of the optical conductivity, 30 ' 31 characteristic size of the observed 
inhomogeneity 9 ' 10 and the doping dependence of the density susceptibility in numerical 27 and 
experimental results, which give ao ~ 0, b ~ 0.7, c ~ 100, F ~ 3X and D s ~ 30X by taking the 
energy unit i(~ 0.4eV) and the lattice constant as the length unit. Basically, all the results 
presented here do not depend on Q within the choice < Q < 0.2-ir. Through the mode 
coupling, the Curie- Weiss type behaviour X( ~ {T + holds in an extended temperature 
region with small Weiss temperature 9 near the marginal critical point (see Fig.2) similarly to 
the spin fluctuation theory. 29 Although the formalism is general, here we adopt carrier doped 
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Fig. 2. The inverse of the density susceptibility peak follows a "Curie- Weiss" temperature dependence 
in an extended region. For the parameter values see the text. 

models with the dispersion of the square lattice E{q) = — 2t(cosq x + cosq y ). The imaginary 
part of the self-energy Im £ is governed by the Curie- Weiss behaviour of xx through eq.(5). 
Although we take ao = 0, the linear temperature dependence is seen in a wide temperature 
region. Therefore, in the marginal-critical region with xx °c X -1 , the resistivity in two- 
dimensional systems becomes nearly proportional to T as p oc ImS oc T. The high-energy 
density fluctuations typically up to around 1 eV also well explains the anomaly observed 
in the eV order of long-tails in the optical conductivity of many compounds near the Mott 
insulator. 2 ' 30,31 

Now we study superconductivity assuming the proximity to the Mott quantum criti- 
cal point. Enhanced density fluctuations mediate the Cooper pairing. The effective inter- 
action between two electrons including that mediated by the density fluctuations is given by 
A(q, iuj n ) = U—U 2 xx(q, i^n)/2 up to the second order in U with xx obtained from (2). To ex- 
tract the role of density fluctuations clearly, we ignore the contribution from spin fluctuations. 
The linearized Eliashberg equation for the superconducting gap A is given as 



for iV-site systems. Here we take a bare Green function ignoring the normal self-energy cor- 
rections to G as the first step. In the calculation, we take the standard Hubbard model only 
with the nearest neighbour transfer t on the square lattice and calculate the eigenvalue A for 
the right-hand side of eq.(6). The self-consistent solution of the linearized Eliashberg equa- 
tion again for the relevant parameter values for the cuprate superconductors has the largest 




k,m 



xA(q - k, i(u n - u m ))A(k, iu m ) 



(6) 
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Fig. 3. The eigenvalue A for the linearized Eliashbcrg equation (7). Among various symmetries, the 
eigenvalue for the d x 2_ y 2 symmetry is the largest and it is plotted as a function of temperature at 
X = 0.1. The transition temperature T sc is given from the temperature of the eigenvalue crossing 
A = 1 (horizontal line). The density susceptibility we used is the same as that in Fig. 2. The local 
interaction in the conventional Hubbard model is taken as U and the nearest-neighbour transfer 
t on the square lattice is taken as the energy unit. 

eigenvalue for the d x 2_ y 2 pairing symmetry and hence the highest superconducting transition 
temperature T sc ( Fig. 3) . Even for small Q, the unconventional pairing comes about with 
the nodal structure due to the attractive part around q = Q coexisting with the repulsion in 
the region far from Q as seen in the form of A. The position of the nodes with the d x 2_ y 2 
symmetry is understood from the largest gap in the (n, 0) and (0, ir) regions stabilized by the 
flat dispersion contributing to the diverging density susceptibility. Even for the bare Green 
function, the d x 2_ y 2 symmetry is stabilized by anisotropic density of states. In Fig. 3, the 
transition temperature has the order of O.Oli to 0.05i, which is the order of a hundred Kelvin 
if we take an appropriate value for the copper oxides (t ~ 0.4eV). 

In summary, we have shown that the quantum Mott criticality causes various effects in- 
cluding strongly momentum dependent structure of quasiparticles, non-Fermi-liquid properties 
and high-temperature superconductivity. We certainly need more elaborate treatment for bet- 
ter understanding on strong-coupling nature of density and doublon-density fluctuations in 
future steps, while all the present results indicate the significance of this proximity. 

The author is grateful to fruitful discussions with K. Kanoda, and D. N. Basov. 
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